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I. INTRODUCTION 



We describe how the thermal counterpart of a vacuum two-point function may be obtained in 
the real time formalism in a simple way by using directly the 2x2 matrices that different elements 
acquire in this formalism. Using this procedure we calculate the analytic (single component) thermal 
amplitude for the pion pole term in the ensemble average of two axial-vector currents to two loops 
in chiral perturbation theory. The general expressions obtained for the effective mass and decay 
constants of the pion are evaluated in the chiral and the nonrelativistic limits. We also investigate 
■ the effect of massive states on these effective parameters. 

(N : 
<: 

The real time thermal field theory is apparently complicated by the fact that all two-point functions in this formalism 
assume the form of 2 x 2 matrices [1] . These matrices, however, have simple structures: If we factor out certain matrices 
depending on the distribution function only, these become diagonal, each with essentially a single independent element 
I , with proper analytic properties. But in actual computations one tends to ignore the matrix structure, starting instead 
' with the so-called physical 11-element, encountering though summation over indices at all interaction vertices in a 
Feynman graph. Such a procedure encounters ill-defined products of components of the matrix propagator, that 
must be combined together to get a well-defined quantity. Further the 11-component does not have a simple analytic 
structure. 

In this work we show that it is both simple and elegant to work with the matrix amplitudes. All one has to do is 
to write out the usual vacuum amplitude. The thermal matrix amplitude is then obtained by replacing its elements 
like the propagator, the self-energy and the vertices by the corresponding matrices. Factorizing these matrices as 
mentioned above, we immediately get the analytic amplitude representing the dynamics of the system in the heat 
bath. 

Here we apply this procedure to calculate the pion pole term in the two-point function of the axial-vector currents 
Qh! to two loops in chiral perturbation theory [2] . This problem was studied earlier by several authors [3-5] , in particular, 
by Toublan [6]. After obtaining the analytic amplitude, we follow him to find the pion pole position and the residue. 
We then find these pole parameters in the chiral limit, in agreement with his results. We also evaluate them in the 
nonrelativistic region and consider the effect of massive states on them. 

In Sec. II we write down the effective chiral Lagrangian to fourth order, needed to obtain all the required vertices. 
In the next Sec. Ill, we obtain the vacuum amplitude from all the Feynman graphs up to two loops contributing 
to the pion pole term of the two-point function. The corresponding thermal amplitude is obtained in Sec. IV, from 
which we derive the effective parameters, namely the pion mass and the decay constants at finite temperature. These 
expressions are evaluated analytically in Sec. V in the high and low temperature limits. We examine the contribution 
of the massive states in Sec. VI. Finally we bring out the main features of our work in Sec. VII. 

Appendix A constitutes an essential part of this work. Reviewing briefly the real time thermal field theory, we 
discuss here at length how the vacuum amplitude for an individual Feynman graph may be converted into its thermal 
counterpart. In Appendix B we write the integrals appearing in the nonfactorizable amplitudes. In the last Appendix 
C we collect the results for the relevant integrals in the high and the low temperature region. 



II. CHIRAL PERTURBATION THEORY 



We consider the QCD Lagrangian for the doublet of light quarks, u and d. In the absence of their masses, it has 
chiral symmetry, being invariant under SU(2)r x SU(2)l. This symmetry is supposed to be broken spontaneously 
to SU(2)v of ordinary isospin, generating the massless pions as the Goldstone bosons. 

In the physical case of non-zero quark masses, chiral symmetry is also broken explicitly to the same isospin subgroup, 
if we neglect the mass difference of u and d quarks. The pions become the pseudo- Goldstone bosons acquiring mass 
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M, given by M 2 = 2m q B to lowest order, where B is related to the quark condensate in vacuum, whose dynamical 
generation leads to the spontaneous symmetry breaking. 

As already stated, we are interested in the pion pole term in the two point function of the axial-vector currents, 

^ = qi^5—q a = 1,2, 3, 

evaluated in chiral perturbation theory, the effective theory of QCD at low energy. Such functions are best calculated 
in the external field method, in which one introduces in the original QCD Lagrangian an external field a£(x) coupled 
to -AjJ(x) as well as a field v%(x) for the vector currents [2]. The global chiral symmetry is then promoted to a local 
one, with appropriate transformation properties of the external fields. 

In the effective theory, the pion fields may be collected in the form of an unitary matrix, 

U(x) = e i +"W T "l F , 

where r a are the Pauli matrices. The constant F may be identified as the pion decay constant in the chiral limit. 
The local symmetry requires us to replace the ordinary derivative by the covariant one, 

D„U = d„U -i{a^U], (2.1) 

where, for our purpose, we retain only the external field a M (x). The two-point function of A^(x) is now obtained 
as the coefficient of the quadratic term in a M (x) in the perturbative evaluation of the generating functional with the 
effective Lagrangian. 

As a non-renormalizable theory, the effective Lagrangian consists of a series of terms with an increasing number of 
derivatives and/or quark mass factors, 

C eff = + £(4) + . . . . . 

The leading term is given by 

£( 2 ) = ^{{D^D»U) + M 2 (U + t/t}} , (2.2) 

where (A) denotes the trace of the 2x2 matrix A. Here the first term is invariant under the chiral transformations. 
The second term represents explicit symmetry breaking due to the quark mass term in the QCD Lagrangian. 
The next, non-leading piece in £ e jj is [2,7] 

£W = h 1 (D^D^U) 2 + -MD^D v U){D^D-U) + \l i M 2 (D^D^U){U+U^) + (l 3 +U)M 4 (U+U^ 2 . (2.3) 
4 4 <3 lo 

It provides counterterms necessary to renormalize the one-loop graphs with vertices from C 1 - 2 ^. Thus the bare coupling 
constants h,--- ,h contain a pole at d = 4 in dimensional regularization. The coefficients of these poles may be 
determined by evaluating all the one-loop graphs. Alternatively, these may be obtained directly by calculating the 
short distance behaviour of the generating functional to one loop [2]. Adopting the notation introduced in this Ref., 
the renormalized coupling constants l\ , ■ ■ ■ ,U are defined by, 

h = 7* f A + T^-Ji) , (2-4) 



with 



32ir 2 



12 1 

7i = o' 72 = o, 73 = -7T, 74 = 2 . 



The pole is contained in A, 



3' u 3' /J 2' 
M d - 4 /ll \ 

A = J^y 2 [ln47r + r (1) + 1] + 0{d 4) ) ■ (2 ' 5) 

Up to the factor 7i/327r 2 , the constants Ij are running coupling constants at the scale M. The M-dependence of Zj 
can be made explicit by relating these to l\, the renormalized coupling constants at any other scale /x, 

'■ = 5^( 7 ' + 1 "f)- < 26 » 

In Sec.V we shall use this equation to show the finitencss of the pion pole parameters in the chiral limit. 

The renormalization of the two- loop graphs would require vertices from the next higher piece, in £ e ff- But we 
do not need it, as we are not interested in the vacuum amplitudes, but only in their temperature dependent parts. 
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III. VACUUM AMPLITUDE 



Here we obtain the pion pole contribution to the vacuum two-point function of the axial- vector current, 

5 a % v (q)=i J ^e^(0|T^(x)4(0)|0) wpole , (3.1) 

by calculating all the Feynman graphs up to two loops with the interaction vertices given by and £( 4 \ These 
graphs are conveniently divided into four groups, shown in Figs. 1-4. 

,,,,,, ,,,-,-Q ,,,,, ,,,,, 
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O ,,„.,,.,,„-,, « ,,,.,,,„, ., : , : , ; .Q O 

(e) (f) 




(h) 



(k) 

FIG. 1. The free amplitude with corrections from one-loop graphs along with counterterm graphs and those two- loop graphs 
that are iterations of the former ones. Vertices of and £* 4 - ) are shown as points and filled circles respectively. Wavy and 
straight lines denote axial current and pion respectively. 
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(d) (e) 
FIG. 2. Factorizable two-loop graphs with vertices from only 
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(d) (e) (f) 

FIG. 3. Further counterterm graphs 
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FIG. 4. Nonfactorizable two-loop graphs 



Let us first derive the familiar one loop results. The vertex correction graphs (b) and (c) of Fig. 1 modify the 
residue of the free amplitude of graph (a) to give 



T { ^ b ' c \q) = q^F 2 {l + 4 V (31 4 - 2J)/3}tA(q) . (3.2) 

M 2 

where r] = is an expansion parameter and A(g) is the free pion propagator, A(q) = i/(q 2 — M 2 + ie). J is a 
divergent one-loop integral, 
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with A given by Eq.(2.5). To include the self-energy graphs, it is convenient to introduce here the well-known Dyson- 
Schwinger equation for the complete propagator A'(q), 

A'(q) = A(q) + A(q)(-iZ(q))A>(q) (3.4) 

where the self-energy part T, of graphs (d) and (e) of Fig. 1 is given by 

E(g) - -2f7(3Z 4 - J)(q 2 - M 2 )/3 + F 2 r/ 2 (4? 3 + J)/2 . (3.5) 

Eq.(3.4) may be solved by iteration, 

A'(q) = A(q) + A(q)(-iZ(q))A(q) + A(q){-iX(q)A(q)} 2 + ■■■ (3.6) 

or in closed form, 

The self-energy correction is now included in Eq.(3.2) by replacing A(q) with A'(q). Thus we get the one- loop result 
for the pion pole, 

= -v? v _^ + ie , 

with 

Ml = M 2 {1 + 2 V (l 3 + J/4)} = M 2 (l - vh/^i: 2 ) , (3.8) 

F 7T =F{1 + r](l 4 - J)} = F(l + riU/16n 2 ) , (3.9) 

on using Eq.(2.4). 

We now include the two-loop graphs. First consider those of Figs. 1, 2 and 3, that are actually products of two 
one-loop parts. The total contribution of all these graphs may be put in the form 

3 

T^ +2+3) (q) = q^q.F 2 ]T{ 7 „ iA(q) + a n M 2 A 2 (q)} + 8F 2 V 2 (h + 2l 2 )(q^ x + q v J» x )q X lA{q) , (3.10) 

n=l 

where we show separately the sums of contributions of all graphs in each of Figs.l, 2 and 3. Thus the sum of graphs 
in Fig.l is given by the n = 1 term 1 . Note here that the two loop graphs from (f) to (1), being iterations of graphs 
from (b) to (e), are automatically included in Eq.(3.7). But we prefer to write them explicitly in Eq.(3.10), getting 

7i = 1 + 2 V {1 4 — J) — 4r/ 2 (Z 4 - J)(3Z 4 - J)/3 , 

g x - r/(4Z 3 + J)/2 - 2r? 2 J(4l 3 + J)/3 , (3.11) 

where we cancel factors as, (q 2 — m 2 )A 2 (q) = iA(q), that can also be justified at finite temperatures. In this way we 
get directly the shifts in the residue and the pole position from graphs of Fig.l as F 2 (ji — 1) and M 2 o\j^\ respectively. 
Next, the graphs of Fig. 2 contain only the vertices of and their sum is given by the n = 2 term, with 

72 = r/ 2 J(8J + 3J')/3 , <7 2 = -?? 2 J(3 J + 2 J')/8 (3.12) 

where we encounter a new divergent integral related to the earlier one as, 

n«) -if**, A> m = -^(m'j) = -^-jt,- <«*> 



: A piece, namely — q ll q 1 y(F 2 /4)ri 2 (J + 4?3) 2 M 4 iA 3 (g), is omitted here, as it is automatically included when we put the n = 1 
term in the form of a simple pole. 
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Lastly, the sum of graphs of Fig. 3 with vertices from and is given by the term n = 3 with 
73 = ?7 2 {(36Zi + 12Z 2 - 25/ 4 ) J + 12hJ' + 12/ 2 }/3 , 

^3 = -r/ 2 {(36Zi + 12Z 2 + I6Z3 - 3U)J + 3l 3 J' + 24(h + 2l 2 )q x q rT J Xa /M 2 } /?> , (3.14) 
together with the remaining term in Eq.(3.10), where we have still another divergent integral, 

J,u(M) = J LJ^L k ,KA(k) . (3.15) 

Actually this term is also proportional to q p q Vl once the integral is evaluated. But in view of its extension to finite 
temperature, we keep it as such. 

Finally we have the amplitude from the nonfactorizable two-loop graphs of Fig. 4, 

T $(>1) = -^U»A( q )T„( q ) + r ii (q)A{q)q v } + J^^A^S^A^) , (3.16) 

where the vertex function r M (g) of graph (a) is 

T » {q) = 1 J ^^ (2 ^ " 3kl » ~ 3 ^)/A(fci)A(fc 2 )A( g -k x - k 2 ) , (3.17) 
and the self-energy function T,(q) of graph (b) is 

£(«) = ~ij ^^( 3M4 + 2/ 2 )A(fc 1 )A(fc 2 )A(g - k, - k 2 ) (3.18) 

with / standing for the function, 

f(q, h,k 2 ) =k{ + kl + Akik 2 + M 2 + 2q ■ (fci + k 2 ) - 2q 2 . 

The order of the factors in Eq.(3.16) is in anticipation of their matrix structures of the thermal amplitudes in the 
next section. 

In Rcf. [6] the vertex and self-energy integrals have been cast in a particularly convenient form using the symmetries 
of the integrands under the interchange of the integration variables. Thus if one defines 

K{q) = JpJ 0|^ A(fcl)A(fc2)A(9 - h ~ k2) (3 - 19) 



they may be written as 



K, v {q) = J ^J ^4 ( rf 27r fc ) 2 4 ^fc 1 ,A(fc 1 )A(fc 2 )A(g - h - k 2 ) , (3.20) 



r„(g) = -2F VM 3J -K) + 9K w qP] , (3-21) 



S( g ) = -F 4 ?7 2 [4(9J 2 - AK){q 2 - M 2 ) + 3M 2 (8J 2 - K) + 72q p q a K pa ] . (3.22) 
Then Eq.(3.16) simplifies to 

- q^F 2 V 2 1 1 J 2 *A(<?) - i(8J 2 -K + 24q"q a K prJ /M 2 )M 2 A 2 (q)^ + 4F 2 V 2 (q^ x + q v K^)q x iA(q) . (3.23) 

The sum of amplitudes (3.10) and (3.23), along with the one from the tree graphs with a single insertion of vertices 
from (not calculated above) would give the complete, renormalized vacuum amplitude. One may then extend the 
one-loop results (3.8, 3.9) for the pion pole parameters to two-loops. Instead, however, we turn to the corresponding 
thermal amplitudes to find the temperature dependence of these parameters. 



5 



IV. THERMAL AMPLITUDE 



The thermal (ensemble averaged) two-point function of the axial-vector current is a 2 x 2 matrix in the real time 
formalism, whose ij-th element is 

(T$)ij = * J d 4 xe^Tr[ P T c Al(Mx)) 4(<^(:c))]/Trp k pole, P = e^ H , (4.1) 

where T c denotes time ordering with respect to the time contour of Fig. 9 in Appendix A. There we discuss at length 
how to obtain the matrix amplitude for an individual graph from its vacuum amplitude. To summarize, all we need 
is to replace the loop integrals (J, J', J^v) encountered in the vacuum amplitudes by (J' 3 , J' @ , J^ v ), where, in effect, 
the vacuum pion propagator is replaced by the 11- or 22- component of the thermal propagator. Further, the elements 
of the vacuum theory, namely, (A, E, T) need be replaced by the matrices (A, X, T) and also A 2 by AtA, where 
the matrices are given in Appendix A. Having obtained the matrix amplitude, we put in the factorized forms for all 
the matrices to get an equation among diagonal matrices, each of whose 11- and 22- elements are identical up to 
complex conjugation and possibly a (— ) sign. Thus we leave behind the matrix structure and work with the (single 
component) analytic amplitude T@ v (q). Almost repeating Eqs. (3.10) and (3.23) we write it as the sum of 

T£< 1+2+3 > = 9 ^F 2 ^{7^A(g) + ^M 2 A 2 (g)} + 8FV(^ + 2Z2)(g M ^A + ?^MA)9 AiA (^ > ( 4 - 2 ) 

n=l 

where 7^ and ct^ are obtained from j n and a n of Eqs. (3. 11), (3.12) and (3.14) after replacing J's by J^'s, and 



l(jP) 2 iA(q) - i{8(J^) 2 -K0 + 24q>>q°KpM 2 }M 2 A 2 (q) 



4F^(q^ x + q v K^iA(q) 



(4.3) 

where, as in Eqs. (3.21) and (3.22) for the vacuum case, we have expressed T^ 3 and T, 13 in terms of K@, K@ v (and J 13 ) 
with 

K " {q) = W>! (t^ (t^ n^OAntfejAn^ - fa - fa) (4.4) 

and 

KM = if, J ff } 4 ( t r fc ) 4 fc ^^Aii(fc 1 )A 11 (fc 2 )A 11 (9 - fa - fa) • (4.5) 

We must point out here that Eqs. (4.4) and (4.5) hold only for the real parts. (The imaginary parts of both sides differ 
by the factor (1 + 2n(|<7o|))~ 1 , as follows from Eqs. (A. 27) and (A. 31).) Since, however, we are interested only in the 
real parts of the pole parameters, our imprecise notation will lead to no error. 

As already stated, the vacuum part of this thermal amplitude is of no interest to us here, beyond the one loop 
results given in Eqs. (3.8) and (3.9). In the /3-dependent part, we have to isolate the finite terms from the divergent 
ones. To this end, we separate the An or . 2 2 into its vacuum and thermal parts in the expressions for J^'s and K@'s. 
In the case of J^'s we have simply, 



where 



K@ splits as 



g, 7 ', 7il „ ) = y|t„ (W) (_L,__^,^) S( e-M* } . (4.6) 



K«{q) = K(q) + K«(q)\ ti , + !( ' + K(g) . (4.7) 



Here the second term gives the /3-dcpcndcnt divergent pieces, all proportional to A. The finite, temperature dependent 
part can be expressed partly in terms of J's defined above and the remainder as certain g-dependent integrals, 
constituting the third and the fourth term respectively. A similar decomposition holds for Kf* v , 

KM = K M + K Ui)Uv +K» + k*»(q) • ( 4 - 8 ) 



6 



All the pieces in Eqs.(4.7) and (4.8) are displayed in Appendix B. 

It is now easy to see that all the /^-dependent divergent pieces cancel out. We then get the complete, renormalized 
thermal amplitude as, 



Te v ( q ) = T^(q) + T, v (q), 

where T' is the vacuum amplitude without the free pole term, which is put in the /3-dependent piece T ^ v 



(4.9) 



T fw (q)=q^F 2 



(1 — 2Jrj + Ar)' 2 )iA(q) + I \jf]— ( B 



K{q) , q\° 



+F 2 {q^S vX {q) + q„S^(q)}qW*Mq) , 
where the tensor is given by 

S^(q) = IJ^ + AK^q) , 
and the (/3-dependent) constants A and B arc built out of J and j' , 



(4.10) 
(4.11) 



A = 7(27 + 0+r(7-A_), n 7:^7. (7 A,). (4,2) 



The three combinations of coupling constants introduced above are, 



1 



12tt 2 



4Z 2 



14 



48tt 2 



6h + 4Z 2 - 9U - 



7 



48tt 2 



6h + 4Z 2 - 6Z 3 - 3/ 4 



55 
12 



(4.13) 



In Eq.(4.10) we have left out regular, non-pole pieces arising out of explicit factors of q 2 . We remove further 
non-pole pieces by expanding the g-dependent functions in qo in the neighbourhood of the pole q 2 = q 2 + M 2 = lo 2 
at fixed q. Thus 



K(q) = K (K "(lu) + (ql - lu 2 )K^(lu) + 



and similarly for S^ u (q). The resulting expression may be put in the form 

U(q).f,Aq) 



TM = - 



q 2 -n 2 (cu)' 



(4.14) 



(4.15) 



where 



1- 



n 2 (to)= q 2 + M 2 { 1 + - Jrj - IB - 



2 ' \ 6 

2 



% U-JV+I A+^r-K (1) (qVS^-J 2 ) A + (W) ( V 



(4.16) 



Note that (q) is not proportional to q^ due to presence of the last term. This is due to the lack of Lorentz invariance 
in a medium, which serves as the preferred frame of reference. Thus, as in nonrelativistic systems [8], we have here 
two different F's, the temporal and the spatial [9] , 



f (q)=q F t (q), 



fi(q)=q l F s (q) 



(4.17) 



One may now find the thermal dispersion curve for the pion and its decay 'constants' at different values of \q\. 
Instead, however, we set q = and find the effective mass and the decay constants as a function of temperature. 
Converting the parameters F and M to the physical values by Eqs.(3.8) and (3.9) and using the result (B.9) of 
Appendix B, we finally get them as 



Ml 



Mi 



11-2 1- 



K 



M 2 {T) = M 2 jl + ^ J -jf(l J+ Y J +I JJ -y +^00 +4^oo 
Ft = F W (1 - #7+** fa + O+ rt + + ^00 + 4ZT 00 - 2M^ 



v 2 ' 2F* 

Mi 



(F*(T) - FZ(T))/F„ = ^ {-12M W C -U + 4(f - K + / J 00 + 4K o)} 



(4.18) 
(4.19) 
(4.20) 
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where l'" is given by 

1 = 19^2 ( 24l l + 1612 - 27l 3 - 2414 - y) ' ( 4 - 21 ) 

and C is the coefficient of the linear term in the expansion of K oi (q) around q = 0, 

C(gb)?i + -". (4.22) 

Note that all the quantities J, K etc. are now functions of M n . These results agree with Ref. [6] except for the 
definition of I [10]. 



V. EVALUATION 



We now need the values of the coupling constants, k, i = 1, • • • , 4. They were already determined in the original 
work [2], but all of them arc not accurate enough [11]. The best values obtained so far follows from matching the 
dispersion theoretic phenomenological representation for nir scattering amplitude to its two-loop evaluation in chiral 
perturbation theory [12], 

Zi=-0.4±0.6, 7 2 = 4.3±0.1, (two loop) 

while the values relevant in the context of one-loop approximation are, 

7i = -1.9± 0.2, 7 2 = 5.25 ±0.04. (one loop) 

The difference in the two sets of values arc attributed to the infrared singularities that can be better dealt with in 
the two-loop matching than in the case of one- loop. The original crude estimate of 73 [2] , 

7 3 = 2.9±2.4,, (one loop) 

has not been improved further. Finally the two- loop estimate of I a [12], 

7 4 = 4.4±0.2, (two loop) 

does not differ much from the original one- loop estimate [2], 

7 4 = 4.3±0.9, (one loop) 

as the infrared singularities are weakly present here. 

It should be noted here that in our two-loop calculation of the pion pole in the axial- vector Green's function at finite 
temperature, it is actually the scattering amplitude in vacuum to one loop that enters its temperature dependent 
part. It is thus appropriate to use the one- loop estimate of the coupling constants in the present context. 

We now evaluate the pion pole parameters in two regions of temperature. First consider the so-called high temper- 
ature limit, T 3> M w . To remain within the domain of the low temperature expansion, this limit is implemented not 
by letting T increase, but instead by holding T fixed and sending M n to zero. The value of M„ is determined by the 
quark masses. Thus the high temperature limit is equivalent to the chiral limit of QCD theory. 

The values of the relevant integrals in the chiral limit are given in the Appendix C. The contributing terms are only 
rfj and the combination, 

r / 2 GJ o + 4i? 00 ) = J^(Z(T)+ C ), (5.1) 

where 

Z(T) = ln^ + l(7 1 + 47 2 ), 
7 

c= -— -ln2 + l-/ 1 + / 2 = 0.30, 
15 

as obtained from (4.13) and (C.l-2). Here Z(T) has actually no logarithmic singularity in the chiral limit, as can be 
seen by shifting the renormalization scale of the coupling constants with Eq.(2.6) from M to any other value fi. 
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We thus get the results for the pion mass and decay constants to two loops at finite temperature in the chiral limit 

as 



M 2 (T) 



Ml 



= 1 + 



= 1 - 



2AF 2 

T 2 
24F 2 

rp2 



12F 2 

T 2 
12F 2 



+ 



Fl(T)-F°{T) 



+ 



36F 4 

36F 4 
72F 4 



In 



A 



M 



Z{T) + c + 



-1 ^ 

72F 4 n T 

f Z(T) + ■ 



where the logarithmic scales are 



A 



M 



1.8 GeV, 



27F 4 V 

T 4 Aa 
27F 4 n T 



A F = A A = 1.6 GeV, 



(5-2) 



(5.3) 



(5.4) 



in agreement with Ref. [6]. Note that A F is associated with F*(T) and not its square, as in this Ref. (Had we 
chosen the two- loop coupling constants, we would get somewhat smaller values for the A's, namely, Am = 1.4 GeV, 
A F = A A = 1.3 GeV.) 

Next consider the low temperature limit, T <C M n . We shall express all quantities in terms of the dimcnsionlcss 
ratio 

T 

times possibly a power of temperature. 

The leading behaviour of all the pole parameters in the low temperature region is given essentially by that of J, as 
seen from Eqs.(C.4-6). Thus we get 



M 2 (T) 



Ml 



\2Fl 24n 2 F 4 J 
M 2 5g 2 M 4 \ - 



F„ 



Fl{T)-F°{T) 



Fl 48tt 2 F 4 ) 
9sM 4 ,. 



\2tt 2 F 4 



(5.5) 
(5.6) 
(5.7) 



where 



27- 3- 9 

gi =h + 2l 2 - —k - + g, 

9 - 3 
g 2 = h + 2l 2 - -r^U + 77T, 



93 



10 
4 

h + 4Z 2 + g, 



10' 



(5.8) 



and J\ T is given by Eq.(C.4). 

It is interesting to note that the two-loop contributions to the effective parameters are always of opposite sign 
compared to that of one-loop, when it contributes. Figs. 5, 6 and 7 show the temperature dependence of these 
parameters. Besides the pion interactions giving these results, massive states are also expected to contribute. But as 
we show below, this contribution turns out to be negligibly small up to T ~ 100 MeV. 



9 



£ 0.5 



25 



50 

T [MeV] 



75 



100 



FIG. 5. Thermal pion mass squared in chiral (continuous curve) and nonrelativistic (dashed one) limits. 
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FIG. 6. 'Temporal' type of thermal pion decay constant in the two limits as in Fig.5. 
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FIG. 7. Difference of temporal and spatial types of pion decay constant in the two limits as in Figs. 5 and 6. 

VI. MASSIVE STATES 

So far we considered only the pionic interactions in obtaining the pion pole parameters. We now examine the 
leading contribution from the state(s) that are massive in the chiral limit. Here the p meson appears to be the most 
important such state contributing through the graphs of Fig. 8. 



(a) w 
FIG. 8. Graphs with p propagator contributing to the pion pole 
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The additional pieces of chiral Lagrangian needed for evaluation of the graphs may be obtained from Ref. [13,14]. 
Their construction follows from noting that the non-Goldstone fields transform only under the unbroken isospin 
subgroup SU(2)v with its group parameters dependent on the Goldstone fields. Accordingly one has to replace the 
ordinary derivative on the p meson field p p (x) by the covariant one, 

V ' p p v = d p p v + [T^p,,]. 

Here the connection T p is given by 

T p = ^u\d p - ia p )u + -u(d p + ia p )u\ u 2 = U, 

where, as in Eq.(2.1), we have set the external vector field to zero. One need also construct appropriate variables from 
D p U and Fj£ L , the right- and the left-handed field strengths built out of the external vector and the axial- vector 
potentials, that transform only under SU(2)y, namely 

u p = iu'D^Uu^ = itj, 

and 

In terms of the above variables, one can write the required pieces of the chirally invariant Lagrangian as [13,14], 

£ mt (p) = —L—(F p { Plu ,f!?) + iG p (p^,u"})), (6.1) 
2v2m p 

where p pv = V ' pPv — X7 u p p . From the decay rates T(p° — > e+e~) and T(p — » 2n), one gets respectively F p = 154 MeV 
and G p = 69 MeV [13]. 

It is now simple to calculate the vacuum amplitudes of graphs in Fig. 8, which may then be converted to their 
thermal counterparts, following discussions in Appendix A. Now the p particles are rather scarce in the medium for 
temperatures up to about 100 MeV [7]. Accordingly we drop the T-dependent term in the p propagator. Further, as 
the momenta involved are small compared to m p , we keep only the leading term in the expansion of the vacuum p 
propagator in inverse powers of m 2 . The amplitudes thus simplified can be expressed through J 13 and J@ v only. 

The self-energy graph contributes a T-dependent amplitude, 

F%q v 24-^ V 2 (J - Joo) (iA(q) + -A 2 (g)J . (6.2) 

Clearly it vanishes in the chiral (high T) limit. Also its leading contribution vanishes in the region where pions move 
nonrelativistically (low T region). In the same way, the vertex graph contributes a T-dependent part, 

F C M 2 

F 2 8 p \ V {-2Joog^ + (q^JuX + J„J„x)q X }iA(q), (6.3) 
m p 

which again vanishes in the chiral limit, but is finite in the low temperature region. Surprisingly, it contributes only 

toF;(t), 



F°(T) 



This contribution, however, is only about 3% of the pionic contribution in Eq.(5.7). 



M ^ F ^ M h\ T (6.4) 



VII. DISCUSSION 



Though the propagators, vertices and self-energies assume the form of 2 x 2 matrices in the real time field theory at 
finite temperature, each of them is essentially given by a single analytic function, as can be seen from an appropriate 
factorization of these matrices. The same is true of the ensemble average of (the T-product of) any two operators. Here 
we show that one can find its thermal 2x2 matrix amplitude directly from the vacuum amplitude and take advantage 
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of this factorization to get the analytic (single component) thermal amplitude. Thus compared to the commonly 
followed practice of considering the 11-clement of the thermal matrix, the use of matrices not only simplifies and frees 
the calculation from ill-defined quantities at intermediate steps, but also yields directly the amplitude with proper 
analytic properties, not possessed by the 11-element. 

In this work we use this matrix method to calculate the thermal pion pole term in the axial- vector two-point function 
in the framework of chiral perturbation theory. From the analytic amplitude we derive the effective mass and the 
decay constants of the pion at finite temperature. These are evaluated in two limits, the chiral and the nonrelativistic. 
The two evaluations agree rather closely up to about T ~ 100 MeV. We also examine the contribution of the massive 
states (p) to the effective parameters and find it to be negligible. 

Finally we compare our work with that of Toublan [6], whom we follow at a number of points. He obtains the 
thermal amplitude in a somewhat intuitive manner, while we formulate rules to write the matrix amplitude, which 
leads immediately to the analytic thermal amplitude. These rules, in effect, justify his way of writing the thermal 
amplitude, as far as its real part is concerned. Our results (4.18-20) for the effective mass and decay constants of the 
pion agree with his, except for the the coefficient V" in the equation for M%(T) [10]. But their chiral limits (5.2-4) 
agree completely with his, as the term with this coefficient does not contribute in this limit. We also find the effective 
parameters in the nonrelativistic limit, in which this term does contribute. In his discussion of the effect of massive 
states on the pion parameters, he considers only the real states in the heat bath and concludes it to be negligible up 
to T ~ 100 MeV. We also consider the effect of the virtual massive states in vacuum and find it too to be negligible. 



APPENDIX A 



The feature of the real time thermal field theory that distinguishes it from the vacuum theory is the time contour 
in their generating functionals. While it is the infinite real line for the vacuum theory, it must be augmented with 
a return path for the thermal theory. One example of such a path that we shall use is shown in Fig. 5. The return 
path may be folded onto the onward path, generating fields with a displaced time argument. Thus compared with 
the vacuum generating functional, 



(0|Texp[ 4 f d 3 x f dt{C mt ((f>) +j(x)(f>(x) +j"(x)A^(x)} in }\0) 

J J —CO 



the thermal one is given by 



Tr 



pTexp[i / d 3 x / dt{C int ((f>i) - C int {<j>2) + ^2 Un(x)</> n (x) + j^{x)A lin {x))} in ] 

J J — oo „_i 2 



(A.I) 



(A.2) 



both written in the interaction representation in terms of the in-ficlds, denoted by the subscript 'in', which we 
shall omit below. The fields 4>i(x) and <f>2{x) have their time arguments on the segments C\ and C2 respectively: 
4>\{x) — (f>(x,t) is the 'physical field' and <p2{x) — (p(x,t — i/3/2) is the 'ghost' field. The — sign before C{4>2) is 
forced upon us by the theory, but the signs before the n = 2 terms with the external fields are at our disposal, which 
we choose to be positive. The pieces C3 and C4 of the complex time contour are of no consequence and have been 
dropped. 



-oo-ip/ 2 



Imt 



+ 00 



Ret 

C 3 



+ oo_£p/ 2 

"-4 
-OO-/3 

FIG. 9. The complex time contour of real time thermal field theory 



The two sets of fields make any thermal two-point function a 2 x 2 matrix. In particular, the free pion propagator 



is 



a«»(x -v)- s« b A(x -v)-( FttWtiiv)) (&iv)<t>iW) \ 



(A.3) 
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where T denotes anti-time ordering, can be evaluated directly in momentum space as, 

_ / A(g) + 2nn(\q \)5(q 2 M 2 ) 2 7 m(|< ?0 |)<% 2 - M 2 )e^ 2 \ 
^ [q) ~ \ 2irn(\q \)5(q 2 - M 2 )e^' 2 A*(q) + 2irn(\q \)5(q 2 - M 2 ) ) ^ AA) 

where n(|<ft)|) = (e /3 l 9 °l — is the pion distribution function and A(q) = ij (q 2 — M 2 + ie) is the pion propagator in 
vacuum. (A boldface letter will always indicate a 2 x 2 matrix.) 

Given the Feynman amplitude for any two-point function in vacuum, it is simple to write the corresponding thermal 
matrix amplitude. The correspondence may be found by comparing the Wick contractions for graphs in the two cases, 
with particular attention to the — sign before the 'ghost' Lagrangian. As an example, consider the graph (d) of Fig. 1, 
for which we show this correspondence in detail. Its vacuum amplitude in coordinate space is obtained from 

~F 2 i J d*z{0\Td^ a {x)d v <ffi(y)Ci nt {<t>{z))\0), (A.5) 

where Ci nt is a piece in 

If-- - - - - M 2 - - 1 

& n t(<t>) - y ■ W -d^-<P- fy# ■ d** -—(</> • ^) 2 } > ( A -6) 

the field <f> denoting the pion iso-vector triplet (<j> , (f> 2 , (f> 3 ). It gives the amplitude in momentum space, 

T™(q) = q»q v J{M)M 2 ||iA(g) + ^M 2 A 2 (g)} , (A.7) 

where J(M) is defined by Eq.(3.3). The single propagator in this expression arises from the cancellation, 

(q 2 -M 2 )A 2 (q) =iA(q). (A.8) 

We now identify the contractions in (A.5) that produce this result. To focus on the contractions, we omit the 
derivatives and isospin indices on the pion fields and write schematically a term of the matrix element (A.5) as 

(O|r0(x)0(2/)0 4 (^)|o) 

~ J{M)A(x - z)A(z - y) 

~J{M)A 2 {q), (A.9) 
in momentum space. To get the corresponding thermal matrix amplitude, consider its ij-th element, 

-F 2 i J d 4 z(Td^Ux)d^ b J (y){£ mt (M^)) - ^nt(M^))})- (A.10) 
Again we write schematically a term of this matrix element and contract its fields as, 

(TfrWMvMiz) -<£(*)}) 

~ JP{M){TU^M{^) ~ <&(*)}) 

~ jP(M){A a (x - z)A lj (z - y) - A t2 (x - z)A 2j (z - y)} 

~/(M)(A(g)rA(g)) 13 , (A.ll) 

in momentum space, where we use in the second line the fact that the contractions of two 4>i's and two 's at the 
same point yield the same result, 

Jfi{M) = w 2 l wr' A{q)llor22 (A - 12) 

and the matrix r is 
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Note that a cancellation similar to Eq.(A.8) for the vacuum case works also here, 

(q 2 - M 2 )A{q)rA{q) = iA(q) . (A.14) 

Comparing the contractions (A. 9) and (A. 11), we see that the thermal matrix amplitude of graph(ld) can be obtained 
from the vacuum amplitude (A. 7) simply by replacing J with J 73 , A with A and A 2 by AtA in it. 

A little reflection on other graphs of Figs. 1-3 will convince us that the thermal amplitudes of all these graphs may 
be obtained from their vacuum amplitudes by the replacements just stated above, together with J' and J M „ by J' 13 
and J@ v respectively, where 

J' f \M) =l J ^ I (A(fc)rA(fc)) llo , 22 
d 



dM 2 



(M 2 J f3 ), (A.15) 



and 



jjLW) = ^4 / k,k„A(k) llor22 (A.16) 
In Eq.(A.15) we use the so-called mass-derivative formula for the matrix propagator [15], 

(M 2 = ^(Ar), (A.17) 

which is the thermal extension of the trivial relation A 2 (q) = idA/dM 2 for the vacuum propagator. 

So long we constructed thermal amplitudes with the matrix propagator only, inserting r explicitly to account for 
the (— ) sign from the 'ghost' Lagrangian. We now introduce two kinds of parts of graphs namely, the self-energy and 
the (two-point) vertices, where it is convenient to include the effect of the associated (— ) sign(s) in their definitions. 
Thus writing the ij-th component of the matrix amplitude of Fig. 4b, again schematically and contracting the fields, 
we get 

~ (A(x - u)S(« - v) A(v - y))ij , (A.18) 
absorbing the (— ) signs in the definition of S, 

E = ( *IL 7" ) > «« - (^) 3 (A.19) 



-S21 S 2 2 



Likewise, for the graph of Fig. 4a we write 



and 



(Tj>Ux){cpt(u) - 4(u)}4>M) ~ (r (1) A)«, r« = ^ j , (A.20) 

(TM*M(«) - <f> A 2 {u)}4>){y)) ~ (AT^j, = ( ^ (A 21) 

The matrix amplitudes are greatly simplified by factoring out matrices involving only the pion distribution function. 
Thus the free propagator given by Eq.(A.4) can be factored as 

A,„ . UW ( a <" A*(q) ) " W ■ " ( ^ ^ ) ■ < A ^> 
Also the full propagator A' and the two-point function T M „ admit similar factorizations, 

A'(?) = t/(g) ( A 'p (,j) A//3 ° ) t/( ? ) , !>(<?) = U(q) ( ^ _r|*( g ) ) ^ ' (A ' 23) 
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as is suggested by evaluation of our graphs. More rigorously, these follow from their spectral representations. 

To derive a similar factorization of the self-energy part £, we look at the Dyson-Schwinger equation for the full 
propagator, 

A' = A + A(-i£)A'. (A.24) 
Inserting the factorizations for A' and A, we infer that £ must have the factorized form [16], 

SM-ir' (*,}«> J.^ )"- 1 - (A-**) 

It immediately follows that 

S 2 2 = -Sn, £21 =£12. (A.26) 
Further we can get the function entirely from £n, 



i?eS /3 = Re "En, ImY? = — l —ImY, ll . 

1 + 2n 



(A.27) 



In the same way the relations, 

T~ -iT (1) A ~ A(-iT (2) ) (A.28) 

give us the factorizations, 

r« = u ( l" % , ) i/- 1 , r( 2 > = t/- 1 ( l" % . ) [/ (A.29) 

We see that and T (2) differ by a — sign in the off-diagonal elements, which is of no consequence to us. We thus 
omit the superscripts to write the relations given by Eq.(A.29) as, 

r 22 = r; 1 , r 2 i = r 12 . (A.30) 

Again the 11 element determines the function T 13 completely, 

ReT 13 = ReT n , IrnT 13 = — - — ImT u (A.31) 

1 + 2n 

APPENDIX B 

The /3-dependent, divergent and finite parts of K°(q) and K& v have been obtained in Ref. [6], which we reproduce 
here for completeness. The divergent parts reside only in terms linear in the distribution function, where we need the 
vacuum integrals, 

L{p) = 1 1 (£^ A(fc) A(p - k) = - 2X ih + R{p) ' 



and 



WAp) = Jpj j^[^k v A(k)A(p - k) 



{(p 2 -6M 2 ). v -4 w 4 



6M 2 



with 



2(24^-M) 2 [2 ° W " 2p2g ^ + {4p » Pu ~ {p2 ~ 4Af2)5 ^ - 4M 2 p^/p 2 }R(p)} (B.l) 
R(P) = -r^2 J dxln (l - A(l - *)/M 2 ) (B.2) 
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Then the terms linear in n are given by 

^4 



KP il)\n = ^J ^^(k 2 - M 2 )n(\k \)L(q k) , 

= - 6Xl - T^ 1+ W J -^S(k 2 -M 2 )n(\k \)R(q-k) (B.3) 



and 



K f ^M\n = j fi J ■0^k,K6(k 2 -M 2 )n(\k o \)L(q-k) + ^ J ^ kfl k v S(k 2 - M 2 )n(\k \)L^(q - k) 

= -^[10J„„ + {5g^ + (4q ll q v - q 2 g^)/M 2 }J] 

! [{(q 2 + M 2 )g^ - 10g M g„}7 - 28M 2 J^} 



2(12ttM) 2 

^4 



+ 6W J ^(k 2 -M 2 )n(\k Q \)R(q-k)- 
[4<fcfc, - 4(g M fc v + g „fc M ) + lO^fe, - ^(q 2 - 3M 2 - 2q ■ k) - AM 2 {q - k)^q - k)J(q - k) 2 ] (B.4) 



Next consider the terms quadratic in n. Introducing the integrals 



Q(p) 

we can write them as 
and 



f d A k 5{k 2 -M 2 )n{k) f s _ f d A k 5{k 2 - M 2 )n(k) 

J (2tt)3 (p-k) 2 -M 2 ' Q ^ p >- J J^ k » ( p -k) 2 -M 2 ' (B ' 5) 

Kf>(q)\ n , =-^J (03^^ - M 2 )n(\k \)Q(q - k) (B.6) 



rfMlni = -^4 J (^3^( fc2 - M 2 )n(|fc |)[{<Z^ + ^k v - 2(q ll k v + g„fc„)}Q(g - k) 

+k lx Q v {q-k) + k v Q lx {q-k)\ (B.7) 

These terms as well as the last terms in Eqs.(B.3) and (B.4) belong to K{q) and K fll/ {q) in the notation of Eqs. (4.7) 
and (4.8). 

Note that Qi and Qo are not independent, but are related by 

Qiip) = ^{2poQo(p) - P 2 Q(P) + J} ■ (B.8) 
Using this relation one can relate and K 00 at the pole, qo = M with q = 0, 

^ = y(^oo-i? +J 2 ) (B.9) 
Being interested in the real parts, we shall not consider terms cubic in n, which are imaginary. 



APPENDIX C 

Here we write the chiral and the nonrelativistic limits of the integrals occurring in expressions for the effective pion 
parameters. The chiral limits were obtained in Ref. [6] 2 , 



2 In Ref. [6] the integral h is evaluated by expressing it in terms of derivatives of Zeta and Gamma functions. But since the 
integral I2 has to be evaluated numerically anyway, we can do the same for I\ also and get identical result. 



1G 



where 



h = 

h = 



15 



dx x 3 In x 

e x - 1 



rp4 



IMF 4 



In- 



M 
T 7 



ln2 + l-/i+/ 2 



(C.2) 



2^ 4 7 e 

is r 00 dxx 3 r 1 

^J e*-lj 



= 0.60, 



da 



1 



(l + a 2 )ln 



1 - a 



1.05 



(C.3) 



The integrals for rj 2 K\ x and ?7 2 §^| x vanish, while those for 'H 2 d Q°° \ x and r/ 2 C| x are finite in the chiral limit. (Actually 
the terms linear and quadratic in n in each of the later two quantities have singular pieces separately in the this limit, 
but they cancel out in their respective sums.) 

Next we calculate the integrals in the low temperature region, where r = T/M v <C 1. Keeping only the leading 
terms, we get 



J 



3 1 



J 



oo 



dK 
dqo 



0{t 





3/2 
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2tt/ 












r 




5/2 e 


-h 



dK. 



00 



9g 



1 Jlr 



16tt 2 



(C.4) 
(C.5) 

(C.6) 
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